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ABSTRACT 

Modifications  were  made  to  Corgnati's  (1971)  steady  state,  symmet- 
rical hurricane  model,  which  was  based  on  a  model  proposed  by  Riehl 
(1963),  and  in  which  bulk  aerodynamic  transfer  formulas  were  used  to 
predict  the  air-sea  interaction  below  a  specified  vortex.   The  modifi- 
cations resulted  in  a  more  flexible  model  in  which  hurricane  intensity 
was  controlled  by  sensible  and  latent  heat  fluxes,  which  were  in  turn 
dependent  on  the  ocean  surface  temperature.  A  boundary  layer  model 
based  on  similarity  theory  was  then  used,  to  specify  air-sea  interaction, 
by  replacing  the  bulk  transfer  formulas  with  Cardone's  (1969)  extension 
of  Blackadar's  (1965)  two-layer  neutral,  baroclinic  boundary  layer 
model.   To  achieve  realistic  results  the  time  scale  was  modified  to 
include  centripetal  acceleration,  which  resulted  in  a  decreasing  height 
of  frictional  influence  with  decreasing  radius.   Varying  Cardone's  for- 
mulation of  roughness  length  and  the  ratio  of  eddy  exchange  coefficients 
for  heat  and  moisture  resulted  in  stress  and  heat  flux  values  in  agree- 
ment with  values  determined  from  hurricane  observations.   Including  the 
complete  boundary  layer  model  resulted  in  greater  sensitivity  to  ocean 
temperature  variations  than  when  bulk  transfer  formulas  were  used  to 
determine  the  air-sea  interaction. 
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An  earlier  thesis  by  McConathy  (1972)  proved  invaluable  in  apply- 
ing and  interpreting  the  results  obtained  from  the  boundary  layer 
model  based  on  similarity  theory. 


I.   INTRODUCTION 

It  is  well  recognized  that  sea  surface  temperature  affects  hurricane 
intensity;  the  hurricane  winds  in  turn,  exert  a  stress  on  the  ocean  sur- 
face and  induce  the  extraction  of  heat  and  moisture.   A  "steady  state" 
hurricane  model  requires  that  a  balance  exists  in  this  energy  exchange 
between  air  and  ocean. 

Riehl  (1963)  described  a  two-layer  model  for  an  axisymmetric  hurri- 
cane in  which  a  balance  existed  between  the  heat  source  and  momentum 
sink  at  the  air-sea  interface.   The  "steady  state"  upper  level  vortex 
was  then  specified  based  on  the  air-sea  interaction  within  the  boundary 
layer. 

Corgnati  (1971)  used  Riehl' s  model  to  construct  a  numerical  model 
which  used  predictive  equations  for  heat  and  moisture  within  the 
boundary  layer,  and  bulk  transfer  formulas  for  fluxes  at  the  air-sea 
interface.   It  was  felt  that  considerable  improvement  in  Corgnati 's 
model  could  be  attained  through  the  use  of  recent  boundary  layer  models 
based  on  Monin-Obuhkov  "similarity  theory";  specifically  through  the 
use  of  Cardone's  (1969)  extension  of  Blackadar's  (1965)  two-layer  neutral 
baroclinic  boundary  layer  model. 

The  purpose  of  this  study  was  to  modify  Corgnati' s  numerical  model 
so  as  to  better  define  the  air-sea  interaction. 

Several  numerical  simulations  of  hurricane  development  have  been 
made  with  complicated  multi-level  models.   These  models  include  some 
interaction  with  the  ocean.   The  scope  of  this  study  was  to  model  only 
the  boundary  layer  beneath  an  upper  level  vortex,  defined  by  Riehl 's 
model,  which  responds  to  the  air-sea  interaction. 
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II.   DISCUSSION  OF  THE  MODEL 

A.   REVIEW  OF  RIEHL'S  MODEL 

Riehl's  two-layer  model  consisted  of  an  inflow  layer  and  an  outflow 

layer  (Figure  1).   Riehl  and  Malkus  (1961)  have  shown  the  conservation 

of  absolute  angular  momentum  is  a  reasonable  first  approximation  in  the 

outflow  layer.   This  constraint  leads  to  the  expression  for  computing 

the  radius  at  which  the  tangential  component  of  the  outflow  goes  from 

cyclonic  to  anticyclonic, 

RQ  =  (2  Vq  R./f)1/2  ,  (1) 

i 

where  R.  is  the  eye  wall  radius  and  V   is  the  tangential  wind  at  R. . 

i 
Riehl  and  Malkus  (1961)  also  demonstrated  the  approximate  conserva- 
tion of  potential  vorticity  in  the  inflow  layer.  Riehl  (1963)  applied 
this  to  the  equation  of  motion  in  polar  coordinates  (see  Figure  2) 
dV    V  V  ,   ..„    ,  dT 

9    J.       Q     r  4-  f\T        ■  1    OP      1      O  fo. 

_  +  __  +  fvr  -  -  -  ^Q  +  -  g-  ,  (2) 

and  integrated  over  the  depth  of  the  inflow  layer  to  get 

RT  =  constant  ,  (3) 

where  T  is  the  tangential  surface  stress.   Using  the  expression 

Te  =  CD  Ps  Ve  /C°S  a  >  (4) 

s 

and  assuming  the  drag  coefficient  (C  )  and  surface  density  (p  )  to  be 

constant,  also  that  the  inflow  angle  is  small  enough  so  that  cosine  cc 

can  be  closely  approximated  by  unit  yields  the  expression 

V^  R1/2  =  constant   .  (5) 
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Figure  1.   Illustrating  hurricane  mass  circulation  V  is  tangential 
wind  component,  9  is  equivalent  potential  temperature,  R  is  radius 


Figure  2.   Illustrating  polar  coordinate  system  used,  V 


wind  component, 
measurement. 


V  is  radial  wind  component. 


R  is  radius 


is  tangential 
9  is  angular 
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Riehl  (1963)  and  Shea  (1972)  have  shown  this  to  be  a  reasonably  good 
approximation  to  the  V  -R  profiles. 

Equation  (2)  can  also  be  solved  for  the  radial  wind  profile, 

where  G  is  the  acceleration  of  gravity,  and  f  is  the  coriolis  parameter. 

Integration  of  the  hydrostatic  equation,  for  typical  values  of 
equivalent  potential  temperature  (9  )  in  a  hurricane,  produces  the 


relationship 


P'  =  -  2.56  0'   ,  (7) 


where  P'  is  the  pressure  departure  from  1005  mb  and  81  is  the  departure 

from  350K.   Riehl  (1963)  also  determined  a  relation  between  the  maximum 

tangential  wind  (V  )  and  the  equivalent  potential  temperature  (6  )  at 

°i  ei 

the  eye  wall  by  equating  V   to  the  height  gradient  of  the  isobaric 

i 
surfaces  and  applying  a  cyclostrophic  relationship  to  obtain 

Ve   (m/sec)  =  14.1(6   -  0  )1/2  ,  (8) 

i  ei    eo 

where  subscript  "o"  denoted  the  value  at  R  . 
r  o 

From  these  equations  it  can  be  seen  that  the  upper  level  vortex  is 
dependent  on  the  increase  in  8  in  the  boundary  layer.   The  problem 
becomes  one  of  modeling  the  boundary  layer  air-sea  interaction  so  that 
the  9  increase  corresponds  to  the  upper  level  vortex. 
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B.   CORGNATI'S  MODEL 

1.   Prediction  Equation 

a.   Atmospheric  Prediction 

In  the  atmosphere  the  predicted  variables  were  heat  and 

moisture  in  the  form  of  potential  temperature  (9  )  and  specific  humid- 

a 

ity  (q  )  respectively.  These  two  parameters,  along  with  pressure  (P), 
determined  9  which  was  the  one  factor  in  determining  the  maximum  wind 
(Equation  8) . 

Partial  derivatives  were  approximated  using  finite  differ- 
encing techniques  to  solve  the  thermodynamic  equation  in  the  following 
form  d9 


EV   9           dV  R0~|      &D0          0  G 
o   a  r     a|  ,      a    ,   _s ,    . 

R^r  +  ^R-j+-^r  +  rA? »  (9) 


oT 

p 

where  cd  is  the  vertical  motion  and  Q  is  sensible  heat  transfer.   The 
first  term  on  the  right  hand  side  was  neglected  due  to  symmetry.   The 
boundary  layer  was  assumed  well  mixed  with  height,  and  therefore,  9 
could  be  calculated  from  the  surface  values  of  air  temperature  and 
pressure.  Moisture  changes  were  determined  in  a  similar  fashion  using 
the  expression 


sr 


pV  q     dV  Rq      W]   Q  G 


where  Q  is  the  latent  heat  transfer, 
e 

Corgnati  (1971)  used  bulk  aerodynamic  transfer  formulas  for 
determining  the  transfer  of  latent  and  sensible  heat  in  the  form 

%  - 1 CP  psVTw  -  Vv  •  (11) 

%  "  f  Lv  psCD(Ow  "  <a)V  '  (12) 

where  C  is  the  specific  heat  of  dry  air,  L  is  the  latent  heat  of 

vaporization  of  water,  p   is  the  atmospheric  density  at  the  sea  surfaces 


T  -  T  is  the  air-sea  temperature  difference,  q  and  q  are 


the 


14 


specific  humidity  of  the  air  and  the  saturated  specific  humidity  at  the 
sea  surface  temperature.   The  surface  wind  (V)  was  approximated  at  70$ 
of  the  wind  at  the  top  of  the  boundary  layer  (Equation  5).   Furthermore, 
Robinson's  (1966)  heat  transfer  coefficient  (2/3  C  )  was  used,  which 
takes  molecular  processes  near  the  surface  into  account. 

Vertical  motion  at  the  top  of  the  boundary  layer  was  com- 
puted from  the  continuity  equation  in  polar  coordinates, 
dV    bRV  v, 

m  +  *sr  +  f?  -°  (13) 

where  the  first  term  is  zero  due  to  symmetry.   Vertical  motion  at  the 
surface  was  assumed  to  be  zero,  and  a  boundary  layer  100  mb  thick  was 
used,  based  on  the  assumption  that  effects  within  the  boundary  layer 
would  not  exceed  this  height. 
b.   Oceanic  Prediction 

Predicted  variables  in  the  oceanic  portion  of  Corgnati's 

model  were  sea  surface  temperature  (T  )  and  the  saturated  specific 

w 

humidity  at  the  sea  surface  temperature.   A  temperature  gradient  of  2C 

per  100  meters  with  a  layer  depth  of  100  meters  described  the  ocean 

structure.   The  amount  of  heat  removed,  through  latent  and  sensible 

heat  transfer,  was  computed  and  the  temperature  of  the  ocean  surface 

reduced  accordingly.   Corgnati  found  that,  since  the  model  was 

stationary,  reductions  in  T  led  to  reduced  air-sea  transfers  and, 

therefore,  reduced  wind  speed.   In  order  to  test  only  the  atmospheric 

portion  of  the  model  the  sea  temperature  was  fixed. 

2.   Initialization  of  Surface  Temperatures  and  Humidities 

The  relationship  between  0   ,  tangential  wind  at  the  eye  wall 

i 
(V  ),  and  pressure  at  the  eye  wall  (P.)  made  it  possible,  by  specifying 
9i  X 


15 


V   ,  to  solve  for  0   and  P.  using  Equations  (7)  and  (8).   The  surface 

ei  ei      L 

pressure  distribution  was  then  determined  by  assuming  a  gradient  wind 

balance. 

Initialization  of  the  temperature  and  moisture  fields,  both 
in  the  boundary  layer  and  at  the  ocean  surface,  were  important  parts 
of  the  model.   The  surface  water  temperature  was  set  at  30C,  then  the 
available  moisture  at  the  ocean  surface  was  calculated  from  the  inte- 
grated form  of  Clapeyron's  equation 

Ln  Oos  "  25  ^f1  -  5-2  Ln  573  •  <14> 

where  es  is  the  saturated  vapor  pressure.   Air  in  contact  with  the 
water  was  assumad  saturated,  with  the  specific  humidity  calculated  from 

.622  e 


P  -  .373e 
s 


(15) 


The  equivalent  potential  temperature  was  computed  by  applying  the 
relationship  between  height  gradient  of  the  isobaric  surfaces  and  the 
tangential  wind  field. 

By  specifying  the  initial  state  from  the  equivalent  potential 
temperature  at  the  eye  wall,  the  effect  of  expansion  cooling  and 
sensible  heat  source  at  the  surface  were  taken  into  account.   Riehl 
(1963)  suggested  that  9  mast  rise  above  350-352K  to  allow  deep  con- 
vection, which  warms  the  column  and  lowers  the  central  pressure  to 
values  which  can  sustain  hurricane  force  winds. 

Specific  humidity  of  the  air  (q  )  was  calculated  assuming  the 
initial  air  temperature  (T  )  equal  to  the  water  temperature,  from 

«.  "  Ta  f   L"  5"  >  (16) 

v      a 

where  T   is  the  pseudo-equivalent  air  temperature. 
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C.   MODIFICATIONS  TO  CORGNATI'S  NUMERICAL  MODEL 

Corgnati  (1971)  used  a  17  by  17  rectangular  grid  with  15  nautical 
mile  spacing.   The  center  of  the  grid  was  the  hurricane  center,  and 
variables  were  computed  along  a  diagonal  and  interpolated  to  the  other 
grid  points.   Due  to  the  symmetry  of  the  model,  the  four  quadrants  had 
identical  results,  and  as  much  information  could  be  derived  through  the 
use  of  a  one-dimensional  model.   Resolution  was  increased  by  using  100 
grid  points  spaced  at  3  kilometer  intervals  along  a  radial.   One  dis- 
advantage of  the  increased  resolution  was  that  smaller  time  steps  had 
to  be  used,  although  this  was  partially  compensated  by  the  reduction  in 
the  number  of  grid  points.   Inherent  instability  in  a  forward -in- time, 
centered-in-space  finite  differencing  technique,  made  it  necessary  to 
change  to  a  leapfrog  technique  with  a  forward  time  step  every  24  time 
steps. 

The  equivalent  potential  temperature  gradient  was  the  determining 
factor  in  predicting  the  maximum  wind.   When  deriving  Equation  (8), 
Riehl  assumed  cyclostrophic  flow,  so  that  the  relationship  was  only  a 
good  approximation  in  the  area  of  maximum  winds.   Assuming  gradient 
winds  led  to  the  relationship 

v2  ,  ,     in3    ae 

^  +  fV  =  -  2'5  x  10   s-S  ,  (17) 

Re        p      dR   '  v  ' 

which  was  used  during  the  initialization  to  determine  the  equivalent 
potential  temperature  field. 

The  interdependence  of  the  equivalent  potential  temperature  field 
and  the  wind  field  makes  the  selection  of  the  eye  wall  and  outer  radii 
important  factors  in  determining  the  wind  field.   The  inner  radius  was 
allowed  to  adjust  towards  the  radius  of  maximum  equivalent  potential 
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temperature,  rather  than  remain  fixed  as  in  Corgnati's  model,  and  the 
outer  radius  was  determined  from  Equation  (1).   The  maximum  tangential 
wind  component  was  then  determined  from  Equation  (17)  using  the  gradient 
of  equivalent  potential  temperature  from  the  eye  wall  to  the  outer 
radius . 

Corgnati  (1971)  found  that  relative  humidity  values  at  the  eye  wall 
region  exceeded  saturated  conditions  after  a  short  period  of  time.   He 
theorized  that  cumulus  clouds,  in  the  form  of  "hot  towers",  with  the 
moist  ascending  air  being  replaced  by  relatively  dry  subsiding  air, 
would  serve  to  remove  this  excessive  moisture.   To  model  the  "hot  towers" 
Corgnati  removed  all  moisture  in  excess  of  a  critical  value  of  95$  rela- 
tive humidity,  and  measured  this  as  convective  precipitation.   A  similar 
problem  was  encountered  at  the  outer  radius,  where  the  computed  value  of 

9  was  found  to  be  significantly  higher  than  the  350K  specified  by  Riehl. 
e 

Since  the  maximum  wind  was  dependent  on  the  radial  gradient  of  9  ,  the 

increased  9  caused  a  decrease  in  the  wind  field.   To  correct  this,  a 
e 

critical  relative  humidity  of  75$  was  used  at  the  outer  radius.   This 
simulated  the  relatively  dry  subsiding  air  outside  a  hurricane.   Several 
profiles  of  critical  relative  humidity  were  experimented  with  and  will 
be  discussed  later. 

Corgnati  (1971)  experimented  with  the  release  of  latent  heat,  within 
the  boundary  layer,  of  the  excess  moisture  removed  through  cumulus  con- 
vection.  Releasing  all  the  latent  heat  within  the  boundary  layer  greatly 
increased  the  hurricane  intensity.   The  bases  of  cumulus  clouds  within 
hurricanes  have  been  observed  to  be  lower  towards  the  eye.   To  model 
this,  latent  heat  release  was  made  an  inverse  function  of  radial  distance. 
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With  these  changes  it  was  now  possible,  through  latent  heat  re- 
lease, relative  humidity  control  and  ocean  temperature  adjustments, 
to  model  steady  state  hurricanes  of  different  intensities. 

D.   CARDONE'S  BOUNDARY  LAYER  MODEL 

The  use  of  bulk  aerodynamic  transfer  formulas  requires  some  know- 
ledge of  heat  and  moisture  exchange  coefficients  as  well  as  winds 
within  the  surface  layer.   Recent  advances  in  the  study  of  atmospheric 
turbulence  have  led  to  a  greater  confidence  in  describing  the  distri- 
bution of  wind  and  temperature  in  the  atmospheric  boundary  layer.   This 
progress  has  arisen  largely  through  the  application  of  the  Monin-Obuhkov 
"similarity  theory",  which  has  been  supported  by  observational  studies. 
Cardone  (1969)  has  extended  Blackadar's  (1965)  two-layer,  neutral, 
baroclinic  boundary  layer  model,  over  fixed  terrain,  to  the  boundary 
layer  over  the  sea.  'McConathy  (1972)  evaluated  Cardone' s  formulation 
of  the  marine  boundary  layer.   The  results  of  McConathy 's  study  were 
used  to  apply  Cardone1 s  specification  of  the  boundary  layer,  instead  of 
the  bulk  aerodynamic  transfer  formulas  in  Corgnati's  model,  and  to 
determine  a  radial  wind  profile  throughout  the  boundary  layer. 

The  Monin-Obuhkov  "similarity  theory"  for  the  constant  stress,  or 
surface  layer,  predicts  that  a  universal  relation  should  exist  for  the 
wind  shear, 

and  the  temperature  gradient, 


bz  z   t 

where  the  friction  velocity  (U^)  is  defined  as 


"  i   ,  (19) 


U*  -  $1/2   ,  (20) 
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and  the  scaling  temperature 

e*  ■  es;   cT  ■  (21) 

*    P  - 
k  is  von  Karman's  constant  of  .35,  z  is  the  height,  H  is  the  total  heat 
flux,  fc  and  <j)  are  functions  of  stability.   Matching  the  integrated  equa- 
tions for  the  Ekman  layer  to  the  integrated  equation  for  the  constant 
stress  layer  at  top  of  the  surface  layer  led  to  the  following  set  of 
equations,  in  terms  of  U^  and  inflow  angle  (a), 

^~=  [2kBQ  sin2a  <>u  (£r)]1/2  ,  (22) 

g 


and 


(23) 


V     LnBQRo  -  f(h/L')   ' 

where  the  integrated  profile  stability  parameter  is  defined  by 

1+i  l+(t2 

t  -  1  -  4>u  -  3Ln^  +  2Ln(-y^)+  2  tan^-jf-LnC-^-Si)  ,  (24) 

and  the  height  of  the  surface  layer  (h)  is  defined  by 

B  V 
h--V  .  (25) 

the  surface  Rossby  number  (Lettau  1959)  by 

Ro-^f   ,  (26) 

o 

and  the  modified  stability  length  (Panofsky  1963)  by 

K.     U^T  ^u/az 

L'  -  r  L  ■  kTWar     •  (27> 

In  these  expressions  V  is  the  surface  geostrophic  wind,  and  B   is 

Blackadar's  constant  of  0.0003.   The  drag  coefficient  (C  )  was  defined 

in  terms  of  the  roughness  parameter  (Z  )  and  stability  parameter  as 

2 


Cz    [L   z/Z  -  Hz/L')J 


(28) 
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The  roughness  parameter  was  formulated  by  Cardone  (1969)  as 

z  =  6'84  S  10   +  4.28  x  10"3  vl   -  4.43  x  10~4  (29) 

°        U*  * 

where  the  second  term  on  the  right  hand  side  represents  Charnock's (1955) 
formulation  for  aero dynamically  rough  flow  under  high  wind  speed  condi- 
tions, ? 

.035UI 

zo=-TJL    •  (30) 

Equations  (22)  and  (23)  were  solved  simultaneously  for  U^/V  and  Oc  from 
the  external  parameters  V  ,  f,  and  air-sea  virtual  temperature  differ- 


E.   APPLICATION  OF  CARDONE* S  BOUNDARY  LAYER  MODEL 

Cardone' s  model  was  applied  to  the  balanced  hurricane  to  determine 
the  fluxes  of  heat  and  moisture  as  well  as  momentum  transfer.   The 
determination  of  a  mean  radial  wind  component,  using  Cardone' s  model, 
assured  continuity  in  determining  advection  within  the  boundary  layer. 

The  tangential  wind,  which  represented  a  gradient  wind  at  the  top 
of  the  boundary  layer,  was  corrected  for  thermal  wind  and  used  as  the 
surface  gradient  wind. 

To  include  the  centripetal  acceleration,  the  coriolis  parameter  was 

replaced  by 

V 
f  =  f  +~    .  (31) 

By  assuming  the  turbulent  transfer  coefficient  for  heat  (K,  )  equals  the 
turbulent  transfer  coefficient  for  momentum  (K  ),  and  knowing  L'from 
the  solution  of  Equations  (22)  and  (23),  the  total  heat  transfer  was 


3 

determined  using      K  U.  C  p  T 

m  *   pa 


^   "         kGH 


(32) 


21 


The  total  heat  (H)  was  then  separated  into  latent  and  sensible  components 
by  assuming  their  ratio  was  determined  by 


Q    c 
-1  =  -£ 


/e    -  e  \ 


The  predictive  Equations  (9)  and  (10)  required  a  radial  wind  for 
advecting  the  heat  and  moisture.   The  assumption  that  heat  and  moisture 
are  homogeneous,  within  the  depth  of  the  boundary  layer,  made  it  possible 
to  use  an  integrated  value  of  radial  wind.  This  was  determined  by 
solving  the  integrated  form  of  Equation  (18) 

o 

at  z  =7   i  =  1,  2,  3,  4,  and  averaging  these  values.   The  mean  radial 
wind  in  the  surface  layer  was  then  determined  by  multiplying  by  the  sine 
of  the  inflow  angle  (cc) .   The  inflow  within  the  spiral  layer  was  deter- 
mined through  integration  of  the  Ekman  solution  through  the  spiral  layer 
to  get         oo 

V  dz=J2   V  since  f— -  [cos(CC4-3£  -Bh)-sin(a+3II  -Bh)],(35) 


Jb. 


where 


■k 

The  two  integrated  values  of  V  from  Equations  (34)  and  (35)  were  then 
used  in  determining  the  vertical  motion  (Equation  13),  and  the  advection 
of  heat  and  moisture  (Equations  9  and  10)  in  the  predictive  portion  of 
the  model. 
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III.   DISCUSSION  OF  RESULTS 

The  numerical  process  encompassed  three  stages.   First  was  an  ini- 
tialization phase  in  which  the  model  was  started  with  a  given  V  ,  T  , 
T  ,  and  R.,  then  the  model  went  through  an  adjustment  process  in  which 
the  moisture  content  in  the  boundary  layer  increased  until  either  a 
balance  was  reached  between  latent  heat  flux  and  moisture  advection,  or 
the  moisture  reached  the  critical  relative  humidity  values.   The  final 
stage  was  a  "steady  state"  in  which  the  time  variations  of  heat  and 
moisture  in  the  boundary  layer  were  small. 

The  following  discussion  deals  with  the  interpretation  of  results 
obtained  when  various  parameters  were  changed,  both  in  the  initial  and 
adjustment  stages  of  the  model.   Due  to  the  complex  interactions  of  the 
model,  a  simplistic  approach  of  establishing  a  base  run  and  then  varying 
one  parameter  at  a  time  was  used.   This  type  of  experimentation  made  the 
effects  of  the  changes  more  obvious. 

A.   CORGNATI'S  MODIFIED  MODEL 

The  modifications  to  Corgnati's  model,  as  discussed  in  the  previous 
section,  were  incorporated  to  achieve  more  flexibility  in  describing 
hurricanes  of  varying  intensity.   Corgnati  discussed  in  detail  the  adjust- 
ment process.   He  found  that  steady  state  was  dependent  on  latent  heat 
release  and  critical  relative  humidity  used  to  simulate  the  cumulus  con- 
vection.  These  parameters  control,  to  a  good  part,  the  heat  and 
moisture  at  the  eye  wall  (R. ) ,  which  in  turn  determines  the  equivalent 

potential  temperature  (6   )  and,  therefore,  the  hurricane  intensity 

i 
(Equation  8) . 
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The  parameters  varied  were  initial  air  temperature,  water  temperature, 
critical  relative  humidity,  and  latent  heat  release  within  the  boundary 
layer.  Figures  3  through  15  illustrate  the  relationships  between  the 
various  parameters,  and  were  obtained  from  the  base  run  with  water  tem- 
perature set  at  30C,  air  temperature  initialized  at  29C,  critical  rela- 
tive humidity  set  at  95$  at  R.,  with  a  linear  profile  to  75$  at  R  ,  and 
no  latent  heat  release  within  the  boundary  layer. 

Air  temperature  (Figure  3)  exhibited  a  nearly  isothermal  (less  than 
,3C)  profile  from  R  to  about  3R.,  with  an  approximately  logarithmic 
drop  from  3R.  to  a  minimum  temperature  at  R. .   The  temperature  drop  cor- 
responds to  the'  very  strong  horizontal  pressure  gradient  (Figure  4) . 
The  adiabatic  cooling  due  to  expansion  would  cause  about  a  5  C  drop, 
but,  as  Riehl  (1954)  pointed  out,  the  blowing  spray  in  the  high  wind 
region  would  greatly  increase  the  effective  surface  area  and,  therefore, 
the  latent  and  sensible  heat  transfer.   This  enhanced  transfer  would  off- 
set the  adiabatic  cooling  and  account  for  the  slight  temperature  decrease 
towards  the  eye  as  observed  by  Deppermann  (1939).   The  bulk  transfer 
formulas  (Equations  11  and  12)  do  not  take  into  account  the  enhanced 
transfer  due  to  the  blowing  spray.   However,  the  increased  wind  speed 
and  air-sea  temperature  difference  were  sufficient  to  increase  the 
sensible  and  latent  heat  fluxes  and  offset  much  of  the  adiabatic  cooling. 

Relative  humidity  (Figure  5)  was  at  the  critical  value  which  indi- 
cates the  steady  state  stage  was  one  in  which  the  flux  and  advection 
of  moisture  was  kept  in  balance  through  the  removal  of  moisture  exceed- 
ing the  critical  relative  humidity.  Figure  15  indicates  the  moisture 
removed  to  keep  the  relative  humidity  at  the  critical  value. 
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Figure  3.   Radial  profile  of  air  temperature  for  the  base  run 
of  Corgnati's  modified  model,  T  =  30C,  T  initialized  at  29C, 
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Figure  4.   Same  as  Figure  3 


except  pressure. 
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Figure  5.   Same  as  Figure  3  except  relative  humidity, 
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Figure  6.   Same  as  Figure  3  except  specific  humidity, 
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Specific  humidity  (Figure  6)  followed  the  relative  humidity  in  the 
isothermal  region.   The  sharp  decrease  in  temperature  reduced  the 
amount  of  moisture,  despite  the  increase  in  critical  relative  humidity, 
and  the  result  was  a  peak  in  specific  humidity  outside  the  eye  wall. 

Equivalent  potential  temperature  (Figure  7)  is  dependent  on  tempera- 
ture, moisture,  and  pressure.   In  the  nearly  isothermal  region  and  with 
small  pressure  gradient,  9   is  governed  by  the  amount  of  moisture.   From 
3R.  to  R.  the  temperature  decrease  has  a  two-fold  effect  in  that  it  also 
reduces  specific  humidity  and  9  becomes  linear.   The  gradient  of  9 

between  R.  and  R  determines  V    (Equation  17),  and  Equation  (5)  de- 

i 
scribes  the  tangential  wind  profile  (Figure  8).   Pressure  and  wind  are 

linked  through  the  gradient  wind  equation,  and  the  sharp  wind  increase 

is  reflected  by  the  strong  pressure  gradient. 

In  a  hurricane  the  winds  rotating  about  the  eye  experience  both 

centripetal  and  coriolis  acceleration  (Figure  9) .   In  the  outer  portion 

of  the  hurricane  coriolis  dominates.   However,  as  V  increases  and 

radius  decreases,  centrifugal  force  becomes  dominant,  supporting  the 

strong  pressure  gradient  in  the  eye  wall  regions,  and  causing  a  decrease 

in  integrated  V  .   The  decrease  in  V  in  the  near  surface  region  requires 

cross  isobaric  flow;  this  results  in  a  frictional  force  opposing  the 

pressure  force.   The  importance  of  these  various  forces  and  their  region 

of  influence  will  become  more  apparent  in  the  discussion  of  results 

obtained  using  Cardone's  boundary  layer  model.   Figure  10  shows  the 

effect  of  cyclostrophic  flow  on  the  radial  wind.   The  strong  vertical 

motion  (Figure  11)  at  the  eye  wall  is  due  to  the  decrease  in  radial  wind 

towards  the  eye  wall  and  the  decreasing  circumference  (Equation  13). 

In  reality  the  vertical  motion  is  due  to  deep  convection  maintained 
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Figure  7.   Same  as  Figure  3  except  equivalent  potential  temperature. 
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Figure  8.   Same  as  Figure  3  except  tangential  component  of  the  wind. 
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Figure  9.   Illustrating  the  forces  acting  in  the  boundary  layer 
of  an  axisymmetric  hurricane. 


Figure  10.   Same  as  Figure  3  except  radial  component  of  the  wind. 
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Figure  11.   Same  as  Figure  3  except  vertical  motion, 
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through  the  buoyant  energy  gain  represented  by  the  9   increase  from  R 
to  R. .   The  convection  is  sufficient  to  cause  the  large  pressure  gradi- 
ent which  is  maintained  by  the  cyclostrophic  balance. 

Latent  heat  flux  (Figure  12)  is  dependent  on  wind  speed  and  the 
difference  in  specific  humidity  of  the  air  and  saturated  air  at  the  sea- 
surface  temperature  (Equation  12) .   The  low  moisture  content  of  the  air 

at  R  causes  a  local  maximum  in  heat  flux.   The  increase  in  moisture 
o 

dominates  the  region  from  R  to  3R.  causing  a  decrease  in  flux  despite 
the  increasing  wind.  Maximum  heat  flux  occurs  at  R.  where  the  decrease 
in  pressure  and  temperature  acts  to  increase  the  specific  humidity 
difference  and  the  wind  reaches  a  maximum  velocity. 

Sensible  heat  flux  (Figure  13),  which  is  dependent  on  wind  speed 
and  air-sea  temperature  difference  (Equation  11),  reaches  a  maximum  at 
R.  and  is  about  an  order  of  magnitude  less  than  the  latent  heat  flux. 

Precipitation  was  simulated  by  two  mechanisms.   One  was  the  large 
scale  moisture  flux  through  the  top  of  the  boundary  layer  (Figure  14), 
which  reflects  the  specific  humidity  and  vertical  motion.  The  second 
mechanism  is  the  removal  of  moisture  which  exceeds  the  critical  relative 
humidity  (Figure  15).   Both  curves  are  accumulated  during  the  first  15 
hours  and,  therefore,  are  influenced  by  the  adjustment  phase.  The  peak 
in  Figure  15  at  R  is  due  to  the  movement  of  R  during  adjustment. 

1 .   Experiments  with  Temperature  Variations 

In  a  series  of  experiments  the  horizontally  uniform  water  tem- 
perature was  varied  from  28  to  31C  while  air  temperature  was 
initialized  at  28C.   The  radial  profiles  of  the  various  parameters  were 
similar  to  those  obtained  from  the  base  run  and  illustrated  in  Figures  3 
to  15.  There  were,  however,  changes  in  magnitude  and  significant  values 
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Figure  12.   Same  as  Figure  3  except  latent  heat  flux. 
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Figure  13.   Same  as  Figure  3  except  sensible  heat  flux. 
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Figure  14.   Same  as  Figure  3  except  total  moisture  flux  through  the 
top  of  the  boundary  layer  for  first  15  hours. 
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Figure  15.   Same  as  Figure  3  except  for  total  moisture  removed  by 
critical  relative  humidity  mechanism  for  first  15  hours. 
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are  presented  in  tabular  form  (Table  1).   Air  temperature  (T  )  at  R 

o 
adjusted  to  within  one  degree  of  the  water  temperature,  and  the  air 

temperature  (T  )  at  R.  was  about  1.5C  cooler  than  T   .   One  of  the 

i  o 

interesting  relationships  was  the  dependence  of  the  model  on  the  differ- 
ence in  9   from  R  to  R. .   Even  when  the  water  temperature  was  fairly 
e       o     x 

cool,  the  increase  in  9  due  to  air-sea  exchange  was  large  enough  to 

sustain  hurricane  force  winds.   Since  9   is  representative  of  the  ambi- 

o 
ent  values  in  a  tropical  atmosphere,  and  values  of  350K  are  necessary 

for  deep  convection,  9   was  limited  to  a  minimum  value  of  350K  for 

o 
determining  V    (Equation  8).   This  limitation  is  reflected  in  the  value 

i 
of  V   for  the  28C  water  temperature  case.   One  of  the  assumptions  of  the 

i 
model  is  that  the  energy  source  for  the  hurricane  is  the  ocean;  there- 
fore, the  hurricane  should  have  more  energy  available  and  increase  in 
intensity  over  warmer  water.   Intensity  did  increase  from  49  m  sec   for 
the  28C  ocean  to  62  m  sec   for  the  31C  ocean. 

2.   Experiments  with  Variations  of  Critical  Relative  Humidity 

Critical  relative  humidity  was  used  to  simulate  decreases  of 
moisture  in  the  boundary  layer  due  to  convective  mechanisms.   To  test 
the  effect  of  the  critical  relative  humidity,  three  profiles  were  speci- 
fied (Figure  16).   Resulting  relative  humidity  profiles  (Figure  17)  show 
that  in  the  convex  case  relative  humidity  did  not  reach  critical  values 
except  in  the  eye  wall  and  outer  regions.   In  the  concave  and  linear 

cases  relative  humidity  reached  critical  values  at  all  radii  between  R 

o 

and  R..   Latent  heat  flux  is  dependent  on  the  difference  in  vapor 
pressure  of  saturated  air  at  sea  surface  temperature  and  ambient  air; 
because  of  this  dependence  the  latent  heat  flux  distribution  was  affected 
by  the  specification  of  critical  relative  humidity  (Figure  18). 


34 


Table  I 

air  temperature  was  initialized  at  28C. 

T      T      T      o      w  e  o 

w      a.     aQ     et     eQ  L 

(°C)    (°C)    (°C)    C°K)    (°K>  (ra/seG)  (mbS> 

28     26.7   27.8   359    347  49  975 

9     27.0   28.5   363    350  7 

10     27.6   29.2   367    353  5y  * 

3?     28.4   29.9   371    356  62  958 
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Figure  16.   Radial  profiles  of  the  three  specified  critical 
relative  humidity  curves. 
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Figure  17.   Radial  profiles  of  relative  humidity  resulting  from 
the  three  different  specifications  of  critical  relative  humidity 


Tw  =  30C,  Ta  initialized  at  29C. 
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Figure  18.   Same  as  Figure  17  except  latent  heat  flux. 
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Figure  19.   Same  as  Figure  17  except  equivalent  potential  temperature, 
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Hurricane  intensity,  which  is  determined  from  8  at  R  and  R. ,  did  not 
J '  e     o      1' 

vary.   The  eye  wall  was  made  to  adjust  towards  the  maximum  9  ;  in  the 
convex  case  maximum  9  occurred  outside  the  eye  wall  and  R.  adjusted  to 

36  km.   In  the  concave  and  linear  cases  maximum  9  occurred  at  or  inside 

e 

the  eye  wall  and  R.  adjusted  to  24  km  (Figure  19). 

It  was  concluded  that  the  actual  form  of  the  critical  relative 
humidity  curve  between  fixed  extrema  was  not  crucial  in  predicting  the 
intensity  of  the  storm;  however,  it  does  affect  the  distribution  of 
heat  exchange. 

3.   Experiments  with  Variations  of  Latent  Heat  Release 

Assuming  dry  adiabatic  expansion  in  the  inflow  would  result  in 
air  temperature  in  the  eye  wall  region  as  much  as  5C  colder  than  ob- 
served temperatures.   Riehl  (1954)  hypothesized  that  latent  heat  release 
or  sensible  heat  flux,  enhanced  by  blowing  spray,  must  offset  the  expan- 
sion cooling.   To  simulate  the  effects  of  this  enhanced  heating,  the 
latent  heat  of  condensation  of  the  moisture  released  to  maintain  relative 
humidity  below  the  critical  value  (Figure  15),  was  added  to  the  boundary 
layer.   Since  the  enhanced  heating  would  be  a  function  of  the  amount  of 
blowing  spray,  and  the  amount  of  blowing  spray  would  be  dependent  on 
wind  speed,  and  wind  speed  is  a  function  of  radius,  the  percentage  of 
heat  released  was  made  a  linear  function  of  radius  with  the  maximum  value 
at  R. ,  and  no  latent  heat  release  at  R  .   Corgnati  experimented  with  re- 
leasing 100$  of  this  latent  heat  at  all  radii  and  found  a  marked  increase 
in  intensity.   Table  II  illustrates  the  effects  of  latent  heat  release  in 

the  modified  model.   The  increase  in  8   ,  while  9   remained  unaffected, 

e.         e 
1  o 

resulted  in  an  increased  intensity. 


38 


Table  II 

Variations  of  significant  values  with  latent 
heat  released  in  the  boundary  layer. 


latent 
release 

heat 

s  at  R. 

T 
a. 

l 

T 
a 
o 

0 
e. 

i 

9 
e 
o 

V     P 

0.     o 

i 

{$) 

(°C) 

(°C) 

(°K) 

(°K)( 

m/sec) (mb) 

0 

27.9 

29.5 

367 

354 

58    963. 

25 

28.3 

29.5 

371 

354 

69    951 

50 

28.5 

29.5 

373 

354 

70    944 

75 

28.7 

29.5 

375 

354 

72    941 

100 

28.6 

29.6 

376 

354 

74    938 
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Although  latent  heat  was  used,  the  heat  gain  could  also  be  pro- 
duced by  an  enhanced  sensible  heat  transfer  due  to  blowing  spray.   In 
all  probability  both  latent  heat  release  and  enhanced  sensible  heat 
transfer  play  some  role  in  offsetting  expansion  cooling  in  the  boundary 
layer.   The  temperature  profile  of  the  model,  without  latent  heat  re- 
lease, was  in  agreement  with  observed  temperature  profiles.   The  model 
was  therefore  run  with  no  latent  heat  release  in  the  boundary  layer. 

B.   CARDONE'S  BOUNDARY  LAYER  MODEL 

Cardone's  two  layer  baroclinic  boundary  layer  model  assumes  hori- 
zontal homogeneity,  stationarity,  K  equal  to  a  constant  in  the  Ekman 
layer  and 

K  =  kU,z/  A  (z/L')  (37) 

m  u 

in  the  surface  layer.   The  model  also  involves  various  empirically 

determined  constants  based  on  data  from  mid-latitudes.   Despite  these 

limitations  it  was  felt  that  a  boundary  layer  based  on  similarity  theory 

would  provide  an  internally  consistent  model  of  an  hurricane  boundary 

layer.   The  inputs  to  Cardone's  model  were  surface  gradient  wind, 

virtual  temperature  difference  between  the  air  in  the  boundary  layer  and 

saturated  air  at  the  sea  surface  temperature,  and  latitude. 

1.   Application  of  Cardone's  Boundary  Layer  Model 

Several  problems  arose  with  the  direct  application  of  Cardone's 

model  in  place  of  the  bulk  transfer  equations  and  the  integrated  inflow 

from  the  tangential  wind  equation.   The  model  hurricane  decayed  rapidly 

during  the  adjustment  phase,  and  heat  and  moisture  fluxes  were  well 

-2    -1 
below  the  3000  to  5000  cal  cm   day   values  determined  by  Whitaker 

(1967)  to  occur  in  the  hurricane  force  wind  region  of  hurricane  Betsy. 

I 
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Investigation  showed  that  radial  wind  (Figure  20),  determined  from  Equa- 
tions (34)  and  (35),  increased  inward  from  3R.  to  R. .   Because  V  must 

11  r 

go  to  zero  in  the  eye,  vertical  motion  at  the  eye  wall  was  on  the  order 
of  1500  mb  hr   .   In  this  region  the  vertical  transport  of  moisture  ex- 
ceeded the  sum  of  the  horizontal  transport  and  evaporation,  which  de- 
creased 8   and  led  to  a  decrease  in  intensity.   Haltiner  and  Martin 
e.  J 

i 

(1957)  have  shown  that  the  balance  of  forces  in  a  vortex  may  be  depend- 
ent on  centrifugal  force.   This  was  included  in  Cardone's  model  by 
replacing  f  by  f  as  formulated  in  Equation  (31).   The  time  scaling 
factor  then  became 

(f  +  V*1)"1  ,  (38) 

which  had  the  effect  of  reducing  the  height  of  frictional  influence  as 
V  R   increased.   In  particular,  the  height  of  the  constant  stress  layer 
was  reduced  from  about  50  meters  at  R  to  about  20  meters  at  R.  (Figure 
21).   To  maintain  consistency  within  the  boundary  layer,  the  centripetal 
acceleration  also  had  to  be  included  in  term  B  (Equation  36)  of  the  Ekman 
spiral.   The  formulation  for  V  within  the  spiral  layer  became 


Vq(z)v  1/2 

\h  u  \2  J  z  sin  a 

where  V  (z)  is  defined  by 


_       /   V  (z)\  1/2 


a  cosfa+^I  f+  R 

L   4  V   u*  k 


(40) 


To  determine  Vr  in  the  spiral  layer  an  iterative  technique  was  employed 
in  determining  Vq(z)  using  heights  (z)  at  10  meter  intervals.   Then  V  (z) 
was  determined  at  each  height  (Figure  22)  using  the  V  (z)  for  that  height, 
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300 


Figure  20.   Radial  profile  of  radial  component  of  the  wind  resulting 
from  direct  application  of  Cardone's  boundary  layer  model. 


Figure  21.   Radial  profile  of  height  of  the  surface  layer,  illustrat- 
ing decreasing  height  of  frictional  influence  due  to  cyclostrophic 
flow. 
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Figure  22.   Vertical  variation  of  radial  component  of  the  wind  in 
the  spiral  layer,  derived  using  iterative  technique. 
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The  integrated  radial  wind  in  the  spiral  layer  was  then  determined  by 
averaging  the  values  of  V  (z)  over  the  spiral  layer.   Using  an  iterative 
technique  proved  to  be  rather  time  consuming.   Because  V  (z)  increased 
rapidly  with  height,  the  assumption  that  V  (z)  was  equal  to  V  at  the  top 
of  the  boundary  layer  resulted  in  integrated  values  of  V  (z)  that  closely 
approximated  the  averaged  values  of  V  (z)  (Table  3).   The  assumption 
that  V  was  constant  with  height  in  the  spiral  layer  was  used  in  deter- 
mining the  integrated  value  of  V   (Figure  23).   Including  centrifugal 
forces  was  now  consistent  within  the  entire  model  and  the  rapid  decay 
described  earlier,  as  well  as  the  extreme  vertical  motion,  was  no  longer 
evident. 

In  a  study  of  water  temperatures  in  the  wake  of  hurricane  Betsy, 

Whitaker  (1967)  estimated  the  heat  flux  in  the  hurricane  force  wind 

-2    -1 
region  to  be  between  .035  and  .058  cal  cm   sec   .   The  application  of 

Cardone's  model  produced  heat  fluxes  with  a  maximum  value  around  .02  cal 

-2    -1 

cm   sec    (Figure  24).   One  of  the  assumptions  in  Cardone's  model  was 

that  the  ratio  of  the  turbulent  transfer  coefficient  for  heat  (IC  )  to 

turbulent  transfer  coefficient  for  momentum  (K  )  was  unity.   Businger  et 

al  (1971)  found  "The  ratio  of  the  eddy  diffusivities  (a.)  for  heat  and 

moisture  is  greater  than  unity.   The  value  at  neutral  is  1.35;  and  shows 

a  marked  increase  with  instability."   Since  heat  transfer,  as  determined 

in  the  model,  is  directly  dependent  on  a,  (Equation  32),  increasing  a, 

n  h 

increased  the  heat  flux.   In  a  hurricane  the  large  quantity  of  blowing 
spray  should  serve  to  enhance  the  heat  transfer.   Furthermore,  the 
velocity  imparted  to  this  blowing  spray  should  serve  to  increase  the 
momentum  transfer.   To  simulate  the  increasing  heat  flux  and  momentum 
transfer,  the  roughness  length  (Z  )  was  increased  by  increasing  the 
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m/sec 
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Table  III 


Comparison  of  integrated  V  values  obtained 
using  iterative  technique,  and  assuming  V  (z) 
constant. 

Radius  (km)   30   60   90   120  150  180  210 

V  (m/sec) 

ilerative    3.14  3.22  3.42  3.49  3.40  3.38  3.20 
technique 

V  (m/sec)        3>22   3^2?   3A7   3^3   3^3   3AQ   3>21 

V  (z)  const 


lio 


R  k. 


300 


Figure  23.   Radial  profile  of  integrated  radial  component  of  the 
wind  assuming  cyclostrophic  flow  and  V  constant  within  the  spiral 
layer.  e 
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Figure  24.   Radial  profile  of  combined  sensible  and  latent  heat  flux 
obtained  using  Cardone's  model,  C  =  .00428,  a,  =  1.0,  T  =  30C, 
and  T  initialized  at  29C.  w 


1 1 /cm2/ 


Figure  25.   Same  as  Figure  24  except  a,  =  2.0,  C  =  .05, 
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2 

the  coefficient  (C_)  of  the  U.  term  in  Cardone's  formulation  of  Z  (Equa- 
2.  *  O 

tion  29).   It  was  found  that  realistic  values  of  heat  flux  were  obtained 
when  C„  was  increased  from  .00428  to  .05,  and  a.  taken  as  2.0  (Figure  25). 

Cardone's  model,  with  cyclostrophic  flow,  C?   increased  to  .05, 
and  a,  increased  to  2,  was  now  used  to  determine  the  sensible  and  latent 
heat  fluxes,  replacing  Equations  (11)  and  (12),  and  the  average  radial 
wind  component,  replacing  Equation  (6). 

A  base  run  was  made  in  which  water  temperature  was  set  at  30C, 
air  temperature  initialized  at  29C,  critical  relative  humidity  set  at 
95$  at  R.  with  a  linear  decrease  to  75$  at  R  ,  and  no  latent  heat  re- 
leased within  the  boundary  layer.   This  base  run  was  compared  to  the 
base  run  illustrated  by  Figures  3  through  15  and  discussed  earlier.   Car- 
done's model  produced  greater  latent  and  sensible  heat  flux  (Figures  12 
and  25)  which  led  to  slightly  warmer  temperatures  at  R.  and  increased 

0  .   As  a  result  of  the  increased  9   ,  V   was  about  5  m  sec   greater. 
e.  e.   o.  6 

i  11 

Table  4  compares  some  of  the  significant  values  of  the  base  runs.   The 
use  of  Cardone's  boundary  layer  model  provided  some  additional  para- 
meters that  were  indicative  of  the  structure  of  the  boundary  layer,  as 
illustrated  in  Figures  26  through  29.   The  modified  stability  length 
(L1),  illustrated  in  Figure  26,  can  be  interpreted  to  be  proportional  to 
the  ratio  of  mechanical  and  buoyant  production  of  turbulent  kinetic 
energy.   The  large  increase  of  L',  with  increasing  V  ,  reflects  an  in- 
crease in  the  mechanical  versus  the  buoyant  production  of  kinetic  energy. 
Surface  stress  values  (Figure  27),  which  increases  rather  dramatically  to 
a  maximum  value  at  R. ,  are  in  close  agreement  with  stress  values  deduced 
by  Miller  (1962)  using  data  from  hurricane  Helene  (Table  5).   The  increase 
in  stress  towards  the  center  is  due  to  the  increase  in  wind  and  the 
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Table  IV 

Comparison  of  significant  values  obtained  from 
two  base  runs. 


Corgnati's 

Cardone's 

parameter 

modified  model 

model 

V    (m/sec) 

58 

62 

T    (°C) 
a. 

27.9 

28.0 

T    (°C) 
a 

29.5 

29.2 

e  °  -  e     (°k) 

e.    e 
1      o 

13 

15 

P   (mbs) 

963 

953 
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Figure  26.   Same  as  Figure  25  except  modified  stability  length. 


Rj   60 

Figure  27.   Same  as  Figure  25  except  surface  stress. 
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Table  V 

Comparison  of  surface  stress  determined 
by  Miller,  and  values  obtained  from  the 
hurricane  model. 


surface  stress 

Radius 

0-20nm 

20-40nm 

40-60i 

dynes /cm 

Miller's 
values 

60.0 

51.0 

27.6 

Radius 

20nm 

30nm 

50nm 

computed 

60 

35 

20 

50 


rougher  surface.   The  increase  of  stress  with  decreasing  radial  distance 
is  reflected  in  the  increase  in  inflow  angle  within  the  constant  stress 
layer  (Figure  28) ;  also  in  the  radial  velocity  at  the  top  of  the  sur- 
face layer  (Figure  29).  These  values  are  comparable  to  inflow  angles 
and  radial  velocities  observed  in  actual  hurricanes.   One  of  the  assump- 
tions made  in  Corgnati's  model  was  that  the  cosine  of  the  inflow  angle 
could  be  closely  approximated  by  unity.  Despite  the  increase  in  the 
inflow  angle,  this  appears  to  be  a  good  approximation. 
2.   Temperature  Experiments  with  Cardone's  Model 

To  further  test  the  use  of  Cardone's  boundary  layer  model,  ex- 
periments were  conducted  in  which  water  temperature  was  varied  from  28 
to  31C.   These  experiments  were  similar  to  the  ones  conducted  with 
Corgnati's  modified  model,  and  discussed  in  section  IIIAl.   The  heat 
flux  distribution  obtained,  using  Cardone's  model,  appeared  to  be  more 
sensitive  to  temperature  variations  than  the  heat  flux  distribution 
obtained  using  the  bulk  aerodynamic  transfer  formulas  (Figure  30).   Due 
to  the  differing  heat  flux,  the  variation  in  the  difference  between 
extrema  values  of  0  was  greater,  using  Cardone's  model,  than  the 
variations  obtained  using  bulk  aerodynamic  transfer  formulas.   The 
result  of  these  greater  variations  was  that,  using  Cardone's  model,  a 

sea  surface  temperature  increase  from  28  to  31C  produced  a  V   increase 

-1  0± 

of  25  m  sec    (Table  6);  while  using  bulk  aerodynamic  transfer  formulas, 

the  same  temperature  increase  produced  only  a  13  m  sec   increase  in 

V  .   It  was  felt  that  the  greater  sensitivity  to  sea  surface  tempera- 

i 
ture,  exhibited  by  Cardone's  boundary  layer  model,  more  closely  simulated 

the  dependence  of  a  hurricane  on  the  oceanic  heat  source. 
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R.    6  0  izu  100  z*»u    p 

Figure  28.   Same  as  Figure  25  except  inflow  angle  in  surface  layer. 


Figure  29.   Same  as  Figure  25  except  radial  component  of  the  wind 
at  top  of  surface  layer. 
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Figure  30.   Radial  profiles  of  combined  latent  and  sensible  heat 

flux,  comparing  values  obtained  from  Cardone's  model  with  T  =  31C 

and  T  =  29C,  and  Corgnati's  modified  model  with  T  =  31C  and 

T  =  W29C,  T  initialized  at  28C.  w 

w       '   a 


Table  VI 

Changes  in  significant  values  due  to  vary- 
ing water  temperature,  air  temperature 
initialized  at  28C. 


(°C)  <°C) 


(°C) 


C°K) 


(°K) 


(m/sec)  (mb) 


28  27.1 

29  27.2 

30  27.8 

31  28.5 


28.0 
28.2 
28.7 
29.3 


358 
363 
368 
373 


348 
349 
351 
353 


45 
58 
66 
70 


979 
960 
952 
944 
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IV.   CONCLUSIONS 

The  results  of  this  study  have  indicated  that  a  marine  boundary 
layer  based  on  similarity  theory  can  be  applied  to  an  hurricane  and 
improve  the  effects  of  air- sea  interaction. 

First,  results  indicate  that  the  time  scaling  factor  in  cyclo- 
s trophic  flow  is  dependent  on  (f  +  V  R  )   ,  rather  than  f  as  in 
Cardone's  model.   The  effect  of  the  centripetal  acceleration  was  to 
reduce  the  height  of  frictional  influence  as  wind  speed  increased  and 
radial  distance  decreased. 

Second,  observed  temperature  fields  indicate  that  enhanced  heat 
flux  is  necessary  to  offset  the  adiabatic  expansion  cooling  within 
the  boundary  layer  of  an  hurricane.   This  enhanced  heat  flux  is  likely 
to  be  due  to  the  large  quantities  of  blowing  spray  observed  in  an 
hurricane.   Consequently,  the  ratio  of  the  coefficients  of  heat  dif- 
fusivity  to  momentum  diffusivity  was  increased  from  unity  to  2. 
Increasing  the  roughness  length,  as  a  function  of  friction  velocity 
squared,  was  used  to  simulate  the  increased  momentum  transfer  and 
resulted  in  surface  stress  values  that  appear  to  be  realistic.  The 
increase  in  roughness  length  was  further  justified  by  resultant  heat 
flux  values  that  were  in  agreement  with  determined  values. 

Third,  the  removal  of  moisture  from  the  boundary  layer,  through 
convective  mechanisms,  appears  to  have  a  strong  influence  on  air-sea 
interaction.   Vertical  advection  by  the  mean  motion  is  not  sufficient 
to  export  the  large  horizontal  flux  of  moisture  within  the  boundary 
layer  and  the  evaporation.   The  removal  of  moisture  in  excess  of  a 
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critical  relative  humidity  value  appears  to  adequately  simulate  "hot 
tower"  convection  and  surrounding  subsidence. 

Finally,  in  the  earlier  model  of  Corgnati  (1971)  it  had  been  neces- 
sary to  specify  a  drag  coefficient  and  extrapolate  a  surface  wind  to 
determine  air-sea  interaction.   However,  the  external  parameters 
required  to  use  similarity  theory  were  surface  gradient  wind,  and  the 
virtual  temperature  difference  between  saturated  air  at  the  sea  surface 
temperature,  and  the  air  in  the  boundary  layer.   The  vertical  structure 
of  the  wind  and  an  implied  drag  coefficient  were  inherent  in  the 
iterative  solution  based  on  similarity  theory. 

The  more  complete  boundary  layer  specification  based  on  similarity 
theory  appears  to  be  justified  in  view  of  the  good  agreement  obtained 
between  predicted  and  observed  or  estimated  values  of  radial  inflow, 
heat  flux,  and  surface  stress.   Furthermore,  the  boundary  layer  model 
calculates  wind  profiles  within  the  surface  layer,  as  well  as  surface 
stress;  and  effects  of  both  stability  and  changes  of  surface  roughness 
on  air-sea  interaction  have  been  included  in  the  hurricane  model. 

One  of  the  assumptions  made  in  this  study  was  that  the  stability 
length  was  defined  by  density  stratification  based  on  virtual  tempera- 
ture.  This  assumption  was  crucial  in  the  determination  of  latent  heat 
flux. 

Further  studies,  in  the  area  of  air-sea  interaction  under  high  wind 
speed  conditions,  are  necessary  to  provide  guidance  in  the  selection  of 
empirical  constants  in  the  similarity  solution.   Specifically,  the 
effects  of  blowing  spray,  in  enhancing  heat  flux  and  momentum  transfer, 
may  be  important  in  maintaining  tropical  cyclones. 


55 


BIBLIOGRAPHY 


1.  Blackadar,  A.  K. ,  1965:   A  simplified  two-layer  model  of  the  baro- 
clinic  neutral  atmosphere  boundary  layer.   Air  Force  Cambridge 
Research  Laboratories  Report  65-531,  pp.  49-65. 

2.  Businger,  J.  P.,  J.  C.  Wyngaard,  U.  Izurai,  and  E.  E.  Bradley, 
1971:   Flux-profile  relationships  in  the  atmospheric  surface 
layer.   J.  Atmos.  Sci.,  28,  181-189. 

3.  Cardone,  V.  J.,  1969:   Specification  of  the  wind  distribution  in 
the  marine  boundary  layer  for  wave  forecasting.   New  York 
University,  School  of  Engineering  and  Science,  Scientific  Report 
GSL-TR69-1,  University  Heights,  New  York. 

4.  Charnock,  H. ,  1955:   Wind  stress  on  a  water  surface.   Quart . 
J.  Roy.  Meteor.  Soc,  81,  639. 

5.  Corgnati,  L.  B.,  1971:   Experiments  with  a  simple  hurricane- 
interacting  ocean  model,  M.S.  Thesis,  Naval  Postgraduate  School, 
Monterey. 

6.  Deppermann,  C.  E.,  1939:   Some  characteristics  of  Philippine 
typhoons.  Manila,  Bureau  of  Printing,  143. 

7.  Haltiner,  G.  J.,  F.  L.  Martin,  1957:   Dynamical  and  Physical 
Meteorology.   New  York,  McGraw-Hill  Book  Co.,  181-183,  234-241. 

8.  Lettau,  H.  H. ,  1959:   Wind  profile,  surface  stress  and  geostrophic 
drag  coefficients  in  the  atmospheric  surface  layer.   Advances  in 
Geophysics,  6,  241-257. 

9.  McConathy,  D.  R. ,  1972:   Specification  of  marine  surface  boundary 
layer  parameters,  M.S.  Thesis,  Naval  Postgraduate  School,  Monterey. 

10.  Miller,  B.  I.,  1962:   On  the  momentum  and  energy  balance  of 
hurricane  Helene  (1958).  National  Hurricane  Research  Project 
Report  No.  53,  10-13. 

11.  Panofsky,  H.  A.,  1963:   Determination  of  stress  from  wind  and 
temperature  measurements.   Quart.  J.  Roy.  Meteor.  Soc,  88,  85-94. 

12.  Riehl,  H. ,  1954:   Tropical  Meteorology.   New  York,  McGraw-Hill 
Book  Co.,  281-356. 

13.  Riehl,  H. ,  and  J.  S.  Malkus,  1961:   Some  aspects  of  hurricane 
"Daisy",  1958.   Tellus  13,  181,  213. 

14.  Riehl,  H.,  1963:   Some  relations  between  wind  and  thermal  structure 
of  steady-state  hurricanes.   J.  Atmos.  Sci.,  20,  276-287. 


56 


15.  Robinson,  G.  D.,  1966:   Another  look  at  some  problems  of  the  air- 
sea  interface.   Quart.  J.  Roy.  Meteor.  Soc,  92,  451-465. 

16.  Shea,  D.  J.,  1972:   The  structure  and  dynamics  of  the  hurricanes 
inner  core  region.   Department  of  Atmospheric  Science,  Colorado 
State  University,  Atmospheric  Science  Paper  No.  182,  Fort  Collins, 
134. 

17.  Whitaker,  W.  D.,  1967:   Quantitative  determination  of  heat  trans- 
fer from  sea  to  air  during  passage  of  hurricane  Betsy,  M.S. 
Thesis,  Graduate  College  of  the  Texas  A&M  University,  65,  28. 


57 


INITIAL  DISTRIBUTION  LIST 

No.  Copies 


1.  Defense  Documentation  Center 
Cameron  Station 
Alexandria,  Virginia  22314 

2.  Library,  Code  0212 
Naval  Postgraduate  School 
Monterey,  California  93940 

3.  Professor  R.  L.  Elsberry,  Code  5lEs 
Department  of  Meteorology 

Naval  Postgraduate  School 
Monterey,  California  93940 

4.  Professor  D.  F.  Leipper 
Department  of  Oceanography 
Naval  Postgraduate  School 
Monterey,  California  93940 

5.  Lieutenant  Commander  N.  A.  S.  Pearson 
Helicopter  Support  Squadron  Seven 
Naval  Air  Station  Imperial  Beach 
Imperial  Beach,  California  92032 

6.  Naval  Weather  Service  Command 

Naval  Weather  Service  Command  Headquarters 
Washington  Navy  Yard 
Washington,  D.  C.   20390 

7.  Lieutenant  Commander  R.  W.  Lyons 
Fleet  Weather  Central 
Norfolk,  Virginia  23511 

8.  Professor  K.  L.  Davidson,  Code  5lDs 
Department  of  Meteorology 

Naval  Postgraduate  School 
Monterey,  California  93940 

9.  Professor  R.  L.  Haney,  Code  5lHy 
Department  of  Meteorology 
Naval  Postgraduate  School 
Monterey,  California  93940 

10.   Professor  R.  J.  Renard,  Code  5lRd 
Department  of  Meteorology 
Naval  Postgraduate  School 
Monterey,  California  93940 


58 


S^cnnty  Classification 

DOCUMENT  CONTROL  DATA  -R&D 

(Security  classification  of  title,    body  of  abstrmct  and  indexing  annotation  must  be  entered  when   the  overall  report  la   classified) 


(Corporate  author) 


Naval  Postgraduate  School 
Monterey,  California  93940 


ia.   REPORT    SECURITY    C  L  A  SSI  F  I  C  A  T  I  Or 

Unclassified 


26.  CROUP 


3   REPORT  TITLE 

Further  Experiments  with  an  Hurricane-Ocean  Interaction  Model 


DESCRIPTIVE   NOTES  (Type  ol  report  and,  Indus  ivo  dales) 

Master's  Thesis;  September  1972 


40RISI  (Firat  name,  middle  initial,  laat  name) 


Nils   A.    S.    Pearson 


6      REPOR  T    O  A  TE 

September  1972 


TOTAL    NO.    OF    Pi 

60 


7b.    NO.    OF    REFS 


17 


•a.    CONTRACT    OR    5R»» 


*.   PROJECT   NO 


iTOH'l    REPORT    NUMSER(S) 


tO.    DISTRIBUTION    STATEMENT 


Approved  for  public  release;  distribution  unlimited, 


I.    SUPPLEMEf 


12.    SPONSORING 


Naval  Postgraduate" School 
Monterey,  California  93940 


13.  ABSTRACT 


Modifications  were  made  to  Corgnati's  (1971)  steady  state,  symmetrical 
hurricane  model,  which  was  based  on  a  model  proposed  by  Riehl  (1963),  and 
in  which  bulk  aerodynamic  transfer  formulas  were  used  to  predict  the  air- 
sea  interaction  below  a  specified  vortex.   The  modifications  resulted  in  a 
more  flexible  model  in  which  hurricane  intensity  was  controlled  by  sensible 
and  latent  heat  fluxes,  which  were  in  turn  dependent  on  the  ocean  surface 
temperature.   A  boundary  layer  model  based  on  similarity  theory  was  then 
used,  to  specify  air-sea  interaction,  by  replacing  the  bulk  transfer  formulas 
with  Cardone' s(1969)  extension  of  Blackadar's  (1965)  two-layer  neutral, 
baroclinic  boundary  layer  model.   To  achieve  realistic  results  the  time  scale 
was  modified  to  include  centripetal  acceleration,  which  resulted  in  a  decreas- 
ing height  of  frictional  influence  with  decreasing  radius.   Varying  Cardone 's 
formulation  of  roughness  length  and  the  ratio  of  eddy  exchange  coefficients 
for  heat  and  moisture  resulted  in  stress  and  heat  flux  values  in  agreement 
with  values  determined  from  hurricane  observations.   Including  the  complete 
boundary  layer  model  resulted  in  greater  sensitivity  to  ocean  temperature 
variations  than  when  bulk  transfer  formulas  were  used  to  determine  the  air- 
sea  interaction. 


DD 


FORM 

I  NOV 

S/N  OJOI -807-681  1 


1473 


(PAGE  1) 


59 


Security  Classification 


Security  Classification 


KEY     WORD* 


Air-sea  interaction 
Hurricane  model 
Boundary  layer 
Latent  heat  release 
Sensible  heat   flux 


DD  All"..  14 73  'back 


S/N    01 Ot  -807-662  I 


60 


Security  Classification 


CP73 


2  1  6  U  0 


Thesis 

138060 

P316 

Pearson 

c.l 

Further  experiments 

with  an  hurricane-ocean 

interaction  model. 

i  ' 

S  E  P  7  3                               216UO 

/ 


Thesis  .;j8060 

P316  Pearson 

c.l  Further  experiments 

with  an  hurricane-ocean 

interaction  model. 


thesP316 

F|™nf,eXPenrTlentS  with  an  hurricane-oc 


3  2768  001  97907  3 

DUDLEY  KNOX  LIBRARY 


